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Abstract 

Knot and link invariants naturally arise from any braided Hopf algebra. We consider the computa- 
tional complexity of the invariants arising from an elementary family of finite-dimensional Hopf algebras: 
quantum doubles of finite groups (denoted D(G), for a group G). These induce a rich family of knot 
invariants and, additionally, arise directly from topological quantum computation. 

Regarding algorithms for these invariants, we develop quantum circuits for the quantum Fourier 
transform over D(G); in general, we show that when one can uniformly and efficiently carry out the 
quantum Fourier transform over the centralizers Z{g) of the elements of G, one can efficiently carry out 
the quantum Fourier transform over D(G). We apply these results to the symmetric groups to yield 
efficient circuits for the quantum Fourier transform over D(5'„). With such a Fourier transform, it is 
straightforward to obtain additive approximation algorithms for the related link invariant. Next, we 
provide efficient quantum circuits for the Clebsch-Gordan transform over D(G) for "fluxon" irreps i.e., 
irreps of D(G) characterized by a conjugacy class of G. For general irreps i.e., those which are associated 
with a conjugacy class of G and an irrep of a centralizer, we present an efficient implementation under 
certain conditions such as when there is an efficient Clebsch-Gordan transform over the centralizers. 
We remark that this also provides a simulation of anyonic models of quantum computation, even in 
circumstances where the group may have size exponential in the size of the circuit. 

As for hardness results, first we note that in contrast to those such as the Jones polynomial, where the 
image of the braid group representations is dense in the unitary group, the image of the representations 
arising from D(G) is finite. This important difference appears to be directly reflected in the complexity of 
these invariants. While additively approximating "dense" invariants is BQP-complete and multiplicatively 
approximating them is #P-complete, we show that certain D(G) invariants (such as D(A„) invariants) are 
BPP-hard to additively approximate, SBP-hard to multiplicatively approximate, and #P-hard to exactly 
evaluate. To show this, we prove that, for groups (such as An) which satisfy two specific properties, the 
probability of success of any randomized computation can be approximated to within any e by the plat 
closure. 



1 Introduction 

Quantum computing using anyons was introduced by Kitacv in his seminal paper [19]. He showed how 
to construct a fault tolerant computer based on anyons where unitary transformations are performed by 
braiding the anyons with each other. The anyons considered in the paper arise from quantum doubles of 
finite groups. Due to the topological nature of braiding, this model is called topological quantum computing. 
Ogburn and Preskill [26] showed how to perform universal quantum computation in this model using simple 
groups. In particular, they showed how to simulate a gate set using the alternating group ^5. Mochon 
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[25] (see also [24]) showed that groups which are solvable but not nilpotcnt can also be used for universal 
quantum computation. 

In a series of papers, Freedman, Kitaev, Larsen and Wang laid the mathematical foundations of topolog- 
ical quantum computing. In [15] and [13], it was shown that quantum computers based on the circuit model 
can simulate topological quantum field theories (TQFT) and hence that topological quantum computers can 
be simulated using the standard circuit model. This, implicitly, gives an algorithm to approximate the Jones 
polynomial of knots and links. Aharonov, Jones and Landau [3] later gave an explicit combinatorial algorithm 
to additively approximate the Jones polynomial. Wocjan and Yard [31] also present an efficient algorithm 
which is based on irreducible representations of Hecke algebras (which give rise to representations of the braid 
group). In the other direction, Freedman, Kitaev and Wang [14] show that topological quantum computers 
based on certain TQFTs can simulate conventional quantum computers. This implies that approximating 
(additively) the Jones polynomial and plat closures of braids is BQP-complete. In [14], BQP-completeness 
was shown when the TQFT was defined at a fixed root of unity k. However, the algorithm in [3] works for 
an asymptotically growing k as well. This gap was closed when Aharonov and Arad [2] showed that even 
when k grows asymptotically, the complexity of approximating the link invariants is BQP-complete. All the 
hardness results are based on one central fact. The representations of the braid group in all these models have 
the property that their image is dense in the unitary group. In [20], Kuperberg showed that this is the case 
in general: whenever the image of the link invariant is dense, additive approximation is BQP-complete and 
multiplicative approximation is SBQP-complete. Using Aaronson's result [1], multiplicative approximations 
are also #P-hard. 

However, very little is known about non-dense invariants such as those arising from quantum doubles 
of finite groups. In this paper, we prove some analogous results for computation with anyons arising from 
quantum doubles of finite groups (denoted D(G')). We provide algorithms and hardness results for approxi- 
mating link invariants in this model. The main difference between the model in this paper and those referred 
to above is that the image of the representations of the braid group in this model is not dense (in fact, it 
is a finite group [12]). This means that none of the techniques developed for dense invariants automatically 
carry over to this model and we develop them in this paper. On the side of algorithms, we first give an 
efficient circuit for the Fourier transform over the regular representation of D(G). We then present quantum 
algorithms to approximate link invariants arising from irreps of D(G). For certain kinds of irreps, namely 
those associated with a conjugacy class of G ("fluxon" irreps), we also give classical randomized algorithms 
to approximate link invariants. Then we give the Clebsch-Gordan transform for fluxon irreps of D(G). For 
general irreps i.e., those characterized by a conjugacy class of G and an irrep of the centralizer, we give an 
efficient Clebsch-Gordan transform under certain conditions. If we can 

1. perform efficient QFT and the Clecbsch-Gordan transform over every Z{g) and Z{g) D Z(h) and, 

2. block diagonalize irreps of centralizers restricted to intersections of centralizers. 

To explain the last condition, note that if p is an irrep of Z{g), then when restricted to Z(g) D Z{h) (a 
subgroup of Z{g)), it breaks up into irreps of the latter group. The last condition, then, says that we must 
be able to perform the transform that block diagonalizes p into blocks of irreps of Z{g) n Z{h). This (rather 
technical) condition can probably be removed. Using these two transforms, we show in the appendix that 
one can simulate anyonic quantum computation efhciently inside even exponentially large irreps of D(G) in 
the circuit model. 

Then turning to hardness results, we show that additive approximations to link invariants for certain 
groups G which satisfy two specific properties are BPP-hard and that multiplicative approximations are 
SBP-hard and exact evaluations are #P-hard. All the algorithms presented in the paper are efficient in 
logdGI) as well as the number of crossings of the link, where |G| is the size of the group. However, our 
hardness results require that the group be of constant size. 

This paper is organized as follows. In Section 2, we introduce quantum doubles of finite groups in terms 
of Hopf algebras. Then in Section 3, we derive their irreducible representations and describe the action 
of R matrices, which are relevant for braiding. We then give an efficient quantum circuit for the Fourier 
transform over D(G). In Section 4, we use this Fourier transform to give quantum algorithms to additively 
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approximate link invariants arising from irreps of D{G). We also give a classical randomized algorithm for 
fluxon irreps. In Section 5, we describe a quantum algorithm for the Clebsch-Gordan transform over D{G). 
Then in Section 6, we show that additive approximations of link invariants are BPP-hard, multiplicative 
approximations are SBP-hard and exact evaluations arc #P-hard. In Section 7, we present our conclusions 
and some open problems. In the appendices, we show how to perform the Clebsch-Gordan transform on 
D(Zp x: Xq) which can be used to perform universal quantum computation. We also show in the appendix 
how to simulate anyonic quantum computation using our transforms, thus using them for a BQP complete 
problem. 

2 Quantum doubles of finite groups 
2.1 Quantum doubles as Hopf algebras 

We recall the notion of Hopf algebra; more complete descriptions can be found in accounts by Kassel [18] and 
Majid [21]. A C- vector space A is a Hopf algebra if it possesses consistent algebra and coalgebra structure 
augmented with an antipode map that yields a natural notion of "inversion." To be more precise, a Hopf 
algebra possesses the following structure: 

Algebra structure A multiplication map fi : A (E) A A and a unit map rj : C A which satisfy 

(Associativity) (g) l)(a 6 (g) c) = (g) /x) (a (g) 6 (g) c) for all a,b,c E A, and 

(Unit) /i(a(g) 1) = (g) a) — a, for all a E A, where 1 is the unique element of A for which ri{c) = cl 
for c e C. 

Coalgebra structure A comultiplication rule A : A ~> A(E) A and a counit e : A ^ C which satisfy 

(Coassociativity) (1 (g) A)A = (A (g 1)A, and 
(Counit) (1 g) e)A = (e (g 1)A = 1. 

Coherence and an antipode The algebra and coalgebra structure are related by the following two axioms: 

(Coherence) The comultiplication map A and the counit e are algebra homomorphisms (where A^A 
is given the natural tensor product algebra structure) and ri{e{a)) — e(a)l. 

(Antipode) An antipode map S : A ^ A satisfying 

H{S (g l)A(a) = ^(1 (g S)A{a) = e(a)l . 

We will denote /i(a g) b) as ab. It follows from the axioms above that the antipode map S : H ^ H is an 
antihomomorphism: 

S{a)S{b) = S{ba). (2.1) 

Furthermore, when H is finite-dimensional (the only case we shall consider), S is invertible as a linear 
operator. 

One perspective on the role played by the coalgebra and antipode structure afforded by a Hopf algebra 
is that it provides ring structure to the family of representations of the underlying algebra. Specifically, 
the comultiplication operator provides the structure of a representation to the tensor product of two rep- 
resentations X and F of A by defining the action of a G A to be that of A (a). The coherence axiom 
A(a)A(6) = A{ab) guarantees that this yields representation structure; the coassociativity axiom yields a 
canonical isomorphism between the representation {X (g) Y) ^ Z and X (E) {Y ig Z). With this notion of 
tensor product of representations, the counit axiom guarantees that the map e : A — > C, a one-dimensional 
representation of A, is a unit under tensor product: X^e = X — eiS)X. Finally, the antipode map gives 
the structure of a representation to Hom(X, C) for any representation X: the action of a G A on an ele- 
ment / S Hom(A, C) is given by {af){x) — f{S{a)x). The fact that this definition defines a representation 
depends on (2.1). 
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A finite group G is naturally associated with two, generally distinct, Hopf algebras. The first is the group 
algebra, the algebra of formal C-linear combinations of group elements denoted C[G'] = {^gCigg \ g S C}. 
The second is the dual of this algebra, the set of maps from G to C. A natural basis for this space is given 
by the delta functions g* , where g*{h) = Sgji- We denote this algebra CG (without brackets). 

Algebra structure on C[G] is given by linearly extending the group multiplication rule; Hopf algebra 
structure on C[G] is defined by adopting the maps 

A{g)=g(E)g, e{g) = 1 and S{g) ^ g-^ . 

Observe that these choices correspond to the familiar notions of tensor product of representations, the trivial 
representation, and representation of dual spaces. It is easy to check that when G is non-abelian, this Hopf 
algebra is non-commutative. Its comultiplication structure, however, is cocommutative, which is to say that 

T-iA(5)T - A(5) , (2.2) 

where T is the linear operator that swaps the two tensor copies A. 

Now consider the dual algebra CG; the algebra structure is given by the rules 

gig2=Sgug2 and 1 = ^ . 

g 

The coalgebra structure is given by the maps 

A(<?*)= J2 92^gl <g*)^Sg,e and S{g*)^{g-'r, 

where e is the identity element of G. For any nonabelian finite group G, CG is a commutative, non- 
cocommutative Hopf algebra. 

With any finite-dimensional Hopf algebra H one may associate a natural dual Hopf algebra H* : each 
map (j> : V W in the definition of H yields a dual map (f>* : Hom(M^, C) Hom(y, C). One can check that 
this process carries the algebra structure of H to the coalgebra structure of H* , the coalgebra structure of H 
to the algebra structure of H*, and that the resulting maps satisfy the axioms. The two algebras described 
above are duals of one another. 

The quantum double D(G) is an algebra defined on the vector space C[G] (E" CG; we write the element 
g ^ h* simply as gh* , letting the superscript on the h remind us that this is an element of the dual algebra. 
The Hopf algebra structure on D(G) is obtained by stitching together the C[G] and CG structures as follows: 

{giK){g2h*2) ^ gig2{hf)*h*2 = S^92 f^^gig2h*2 and 1 = e^h* 

h 

where = y^^xy denotes the conjugation of x by y and 

A(gr)- ^ gh*2®ghl, e{gh*)^6h., and S{gh*) ^ {g-'){gh~'g''r . (2.3) 

hih2—h 

Remarks. If is a Hopf algebra with comultiplication rule A, one obtains the coopposite of this algebra, 
denoted H^°'^ ^ by reversing the comultiplication rule: specifically, the comultiplication structure of iJ'=°P is 
given by r o A, where T: H®H^H®H\s the exchange operator (that linearly extends the rule T : 
a(E)/3 I— )■ /3(E) a). All other algebraic structure is inherited from H . In fact, it is this reversed comultiplication 
structure of CG'^°p featured in the definition (2.3) of the quantum double construction. D(G) is generated, as 
an algebra, by the elements {^f^gh* \ g £ G}, an embedded copy of C[G], and {\h* \ h e G}, an embedded 
copy of CG'^"?. 

If G is non-commutative, then the quantum double D(G) is a non-commutative, non-cocommutative 
Hopf algebra. The construction is a rough analogue of the semidirect product construction of two groups: 
specifically, observe that one can interpret the multiplication rule 

{9ihi){g2h*2) = gig2{hf)*h*2 = S^i^ ,^^gig2h*2 
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as a consequence of the familiar "commutation relation" hg = g{h^) where, as above, we adopt the notation 
= ghg~^. Physically, one can think of the elements of the group as creation operators and the elements of 
the dual as annihilation operators. The multiplication rule, then, expresses a natural commutation relation 
typically more exotic than the usual one for bosons {aa^ — a^a) or fermions {aa^ — —a^a). 

Quasi-triangularity; braiding Quantum doubles of finite groups possess a further property of interest: 
they are quasi-triangular (or braided) . Specifically, there is an invertible element R oi A(>^ A that imparts a 
"near cocommutativity" property on A, 

RA{a)R-^ = TA{a) , (2.4) 

and satisfies the braiding relations on D{G)'^^ = D(G) <g) D(G) D{G): 

{A<»l){R) = Ri3R23, (1® A)(i?) =i?i3i?i2, (2.5) 

where R12 (and R23, respectively) denotes the element R acting on the first two (last two, respectively) 
tensor copies. (As above, T:a(8>/3H'/3(8)ais the exchange operator.) 

The element R is traditionally called an "i?-matrix," and immediately yields a solution to the Yang-Baxter 
equations: if we define si = TR12 (E" / (and likewise for S2 and S3), then it can also be shown from (2.5) that 

S1S2S1 = S2S1S2 . (2.6) 

In the algebra 0(0)®", we generalize this notation by analogously defining s^, for 1 < i < rt, as / (g) TR (g) /, 
the operator TR acting on the zth and i + 1st factors of D{G). 

A generalization of the quantum double construction described above for finite groups can be applied to 
generic Hopf algebras (and their dual) to yield a quasi-triangular Hopf algebra with an explicit expression 
for the R matrices. For the finite group double D(G), the R matrix is given by 

R = J29^9* e D(G)® D(G), 
9 

where we embed g E G into D(G) as J2h9^* ^^"^ 9* ^9* ■ 
2.2 The braid group 

Recall that the braid group on n strands, i?„, is an infinite, discrete group generated by the elements 
{cti, . . . , dn-i} with the following relations: 

aiffj = cTjO-j for \i~ j\ > 2 , 

ai(7i+iai = (Ti+icr,;cri+i for i G [l,n - 1] . (2.7) 

With the natural topological realization of the braid group as equivalence classes of braids under ambi- 
ent isotopy, the generator ai can be naturally associated with the topological braid of Figure la and the 
relations (2.7) correspond to the topological equivalence of Figure lb. 

Representations of the braid group via the quantum double As the elements Si of (2.6) above 
satisfy the relations (2.7), it follows that if p : D(G) — > GL(y) is a representation of D(G) we can obtain a 
representation t = p*^" of i?„ on V^"- by letting T(cri) = p'^"(sj). 

3 The Quantum Fourier transform over D{G) 
3.1 Background in representation theory 

We assume that the reader is familiar with the representation theory of finite groups and merely set down 
notation below, following Serre [29]. 
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(a) The topological realization of the element Oi, in- (b) The topological interpretation of the Yang-Baxter 
volving strands i and i relations: aiai+iai — ai+iaiai+i. 



Figure 1: Topological realization of the braid group. 



Let G be a finite group. A representation p of G is a homomorphism p : G — > U(F), where is a finite- 
dimensional complex Hilbert space and U(V^) denotes the unitary group on V. In this unitary case, it is 
convenient to work with a G-invariant inner product, one for which (u, v) = {p{g)u^ p{g)v). By symmetrizing 
an arbitrary inner product (•, •) on V , one can define a new, G-invariant inner product by the rule 

{u,v)' = ■^^^{p{g)u, p{g)v) ; 
' ' see 

thus we work under this assumption without loss of generality. 

To be concise, we will simply say that V is a representation of G (omitting reference to the map p) and 
let gv denote p{g)v, the linear action of an element g € G on a vector v G V. Two representations V and V 
are isomorphic if there is a linear isomorphism (f> : V ^ V for which g(j){v) = 4>{gv) for all g and v. We say 
that a subspace of a representation V is invariant if gw G W for all g G G and w € W. Such an invariant 
subspace is a representation itself under the restricted action. Of course, both {0} and V are invariant; if 
these are the only invariant subspaces, V is irreducible. When V is not irreducible, there is a nontrivial 
invariant subspace W C V and it is easy to check that W'^ is likewise invariant. This realizes V as the direct 
sum of two, orthogonal invariant subspaces, each a representation of G. When a representation p : G — ?> U(l^) 
can be decomposed in this way so that V = Wi (B W2, we write p — (Ji (B cr2, where at : G ^ U(lVi) is the 
restricted action on Wi. This process allows one to decompose any finite-dimensional representation into 
an orthogonal direct sum of irreducible representations. While this decomposition is not, in general, unique, 
the number of appearances of representations of each isomorphism class is determined uniquely. 

For a finite group G, the vector space C[G] can be naturally given the structure of a representation by 
linearly extending the left multiplication map g : h 1-^ gh; this is a representation of degree |G| called the 
regular representation which we will alternately denote Regg. Likewise, the one-dimensional vector space 
C can be given the structure of a representation by linearly extending the rule 17 : 1 i~> 1; this is the trivial 
representation. The regular representation C[G] has a remarkable, canonical decomposition into irreducible 
representations of G: each irreducible representation of G appears with multiplicity equal to its dimension. 
We write 

C[G]=0p®^% (3.1) 
p 

where this sum extended over all (isomorphism classes of) irreducible representations p : G — > U(V^) and dp 
denotes dim Vp. 

A (unitary) linear operator on C[G] that carries the basis {g \ g E G} to a basis consistent with the 
decomposition (3.1) is called a Fourier transform (over G). When convenient, we adopt the Dirac vector 
notation to avoid lengthy subscripts and write target basis for ^pP®'*'' in the form \p,i,j), where p is a 
representation of G, and i,j £ {1, . . . , dp}. Here we treat j as indexing one of the dp copies of p in C[G]; 
thus, for each p and j, the space spanned by {|p, i, j) \ i G {1, . . . , dp}} is an irreducible subspace isomorphic 
to p. 



6 



Given a subgroup H oi G and a representation {p, V) of H, one can construct a representation of G 
from p by a process called induction which plays an important role in the representation theory of D(G). In 
this subsection, we describe the construction in general. As a G module, the induced representation can be 
neatly expressed as a tensor product (over the noncommutative ring C[H]) 

C[G] (g>ciH] V, 

where ®c[h] indicates that this tensor product is C[i/]-linear, i.e., gh(E)v = g ^ p{h)v, for h £ H and v £ V. 
In order to give an explicit construction in terms of a basis, let us pick a transversal T = {ti, t2, ■ ■ ■ ,tk} for 
H in G (so that k — \G\/\H\) and fix a basis {\vi), . . . , \vd)} for the representation space V. The induced 
representation, which we denote p or simply f*^ p when H is understood, acts on the vector space C[T]^V, 
where C[T] = I ij € C} is the (Hilbert) space of formal C-linear combinations of elements of T 

with the inner product {U \ tj) = 5ij. (The tensor product in this case is the conventional tensor product 
over C.) In this basis, (t/| p){g) is given by linearly extending the rule 

g : \t) (g) \v) ^ \t') ® h\v) (3.2) 

where t' £ T and h £ H are the unique elements for which gt — t'h. 

There are two facts about induced representations which we use later. 

1. Induction is transitive: ii H < K < G then 

[tl(tgp)]=tgp. (3.3) 



2. The regular representation of G is an induced representation from the trivial representation of the 
trivial subgroup {e}. 

RegG =tf,} 1 . (3.4) 

Combining the above two facts, we can write 

RegG =tg (tfe} 1) -tg Reg^. (3.5) 

3.2 Irreducible representations and Fourier analysis over 0(6*) 

A representation of a Hopf algebra H is a representation of its algebra structure. Specifically, a representation 
(p, V) of the Hopf algebra is a homomorphism 

p:H^ GL{V) . 

As above, the algebra D(G) can be given the structure of a representation using left-multiplication. The 
irreducible representations of D(G) have been worked out in Dijkgraaf et al. [11] in the physics literature and 
by Gould [16] in the mathematics literature. In the article [16], Gould established that the algebra D(G) is 
semi-simple (when G is a finite group), described the irreducible representations of D(G), and developed an 
associated character theory. In this section, we briefiy describe his construction of the irreducible representa- 
tions of D(G), indicate how the regular representation can be broken down into irreducible representations, 
and proceed to construct quantum circuits for the associated quantum Fourier transform. 

To simplify the notation in what follows, we adopt Dirac notation for the algebras C[G] and CG, treating 
them as spanned by the orthonormal bases {\g) \ g £ G} and {\h*) \ h £ G}, respectively. Writing \gh*) as 
shorthand for the tensor product \g, h*), the action of D(G) on the regular representation, spanned by the 
\gh*), is 

g'(h'y\g,h*)=5^h,y,h\g'g,h*). (3.6) 

We see that g'{h')* takes \g,h*) to zero unless h and h' are conjugates (under the action of g); in the case 
where g'{h')* does not annihilate \g,h*), it acts by the left G action in the first coordinate. It follows 
immediately that each of the subspaces 

Vh^span{{\g,h*),g£G}) 
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is invariant under the left action of D(G); this provides an initial orthogonal decomposition D{G) = 0^ Vh 
into invariant subspaces. Observe, also, that the action of X^^eg 9^* on Vh is precisely left multiplication 
by g on the first component; this gives Vh the structure of a G-representation isomorphic to C[G]. As 
the action of D(G) on Vh is at least as rich as the action of C[G] on Vh, the irreducible subspaces of 
Vh under the D(G) action are direct sums of irreducible subspaces under the G action.^ Remarkably, we 
will see that these representations arise directly as induced representations from the centralizer subgroup 
Z{h) ^{gGG\gh^ hg}. 

Fixing an h G G, let Th be a transversal for Z(h) in G; then we may write each basis vector of Vh 
uniquely: 

\g,h*) = \tz,h*), 

where t € Th and z G Z{h). In particular, we have the straightforward isometric isomorphism C[G] = 
C[T;i] (g) C[Z{h)] given by the map \g) \t) (g) \z). With this isomorphism, we shall work with the basis 

{\t,z,h*)} (3.7) 

(for Vh). Recall from the discussion of induction, above, that C[G] = tz{h) ^^Sz{h)- Indeed, this \t,z,h*) 
basis is precisely the basis described in the explicit formulation of induction given in (3.2). 

It is easy to check that induction commutes with decomposition: when iJ is a subgroup of G and p and 
T are two iJ-representations we have a natural isomorphism between the two G-representations f/f {p © ''") 
and til P ffi til According to (3.1), the group algebra C[Z(/i)] decomposes into the orthogonal sum 
0^p®'^'', in which each irreducible representation p of Z{h) appears dp — dimp times. We conclude that 
C[G] = 0p tfft t'^is sum extended over all irreducible representations of Z{h). We remark that when 
p is an irreducible representation of Z{h), the representation t^(^h) P "ciot, in general, irreducible as a G 
representation. We shall see, however, that the natural D(G) action on '\'z[h) P irreducible; in fact, all 
irreducible D(G) representations arise in this way. Carrying out the Fourier transform on the component 
corresponding to Z(h), we finally arrive at the basis |i, p, h*)\ for convenience, we reorder the components 
and work with the basis 

{\h\p,j,t,i)]. (3.8) 

For a given p, h, and j, the vectors \h* , p, j,t,i) span an induced representation of p to G; we denote this 
space as Vh^pj. The action of G on these vectors is that of t'z{h) P'^ ^'^ particular, G leaves this space invariant. 
Though the space is not irreducible-in general-under the G action, it is under the richer action afforded by 
D(G). In light of (3.6), the action of any {h')* preserves this space; specifically, 

• if h' is not a conjugate of /i, it annihilates Vh and hence Vh,p,j] 

• if h' is conjugate to h, the action of h* on the space Vh is the projection onto the span of 

{\g,h*)\h' ^h<^}. 

Note, however, that /i^^ — for any z € Z(h) and hence {g\h' = g''} is a left coset of Z{h). Evidently, 
there is an element t' e Th for which h*\t, z,h*) = (5^ |i, z, ft,*) (in the basis of (3.7) above) and we 
conclude that on the space Vh,p,j the action of h* projects onto the span of the vectors 

{\h,p,j,t',i) \ ie {!,..., dp}}. 

We let p J denote this subspace. 

Note that Vh^pj = ^teT^ pj' orthogonal direct sum. To see that this space is irreducible under the 
D(G) action, consider a D(G)-invariant subspace W C Vh^pj. As described above, the D(G) action contains 
the orthogonal projection operators onto each Vj^ p^; it follows that 

= Wt where each Wt = W n V^pj . 

^The map g i— > gh* is a one-to-one algebra homomorphism from C[G] into D(G); thus any D(G) representation can be 
given the structure of a G representation. 
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Now we turn our attention to the G action. For concreteness, we assume that 1 £ is the representative 
for the coset Z{h); then note that the subspace Vh p j is invariant and, moreover, irreducible under the Z(h) 
action as it is isomorphic to the irreducible Z(/i)-representation p. It follows that Wi is either or p j- As 
Vfi.pj has the structure of an induced representation (it is isomorphic to Ind^(;j-|p), the G action is transitive 
on the spaces V/ipy. in particular, for any ti,t2 G Th, there is an element g £ G so that g{Tl^pj) C T^^pj 
and, hence, g{Wti) C giWt^)- We conclude that dim Wt^ = diml^^t^ foi' ^^^h ti,t2 and, finally, that W = 
or Vh,p,j, as desired. 

This yields an orthogonal decomposition 

D(G)=e e e v,,, 

into irreducible subspaces. Here the notation Z{h) denotes the collection of irreducible representations of Z{h) 
(upto isomorphism). As D(G) is semisimple [16], all irreducible representations appear in this decomposition. 
This decomposition will be sufficient to describe the quantum Fourier transform over 0(0). 

To complete the representation-theoretic picture, however, we discuss a few more details. It is not difficult 
to show that the D(G) representations V^.pj and Vh' ,p' j' are isomorphic whenever the pair (h, p) and {h' , p') 
are conjugate in the sense that there is a group element g so that h' ~ and the Z{h') representation 
(p')f : z p'{z^) is isomorphic to p. 

In order to describe an explicit action of D(G) on its constituent irreducible representations, note that 
for a transversal we can pick elements kg such that {kg)^^gkg — h, where g is a conjugate of h. The the 
action of D(G) can then be written as 

gh*\kg',v) ^ 5h,g'\kgg>g-i,p{k~^,g_^g'kg,)v) . 

For notational simplicity, we denote the transversal element kg by g in the first register and write the above 
equation as 

gh*\g',v) = 5hMg-\p{k-J,g-,g'kg,)v) . (3.9) 

Finally, a remark about conjugate (contragredient) representations over D(G). If p : D(G) — > GL(y) 
is a representation of D(G), its conjugate representation p* : D(G) — > GL(y*) is defined by the rule 
p*{a)f : V I— >■ f{p{Sa)v). It follows that the conjugate of the irreducible representation liidz(h)P is the 
representation lndz{h-^)P* ■ Self-dual representations will play a special rule in our applications to knot 
theory; see Section 4.1.2. 

3.3 The quantum Fourier transform over D{G) 

The quantum Fourier transform (QFT) over D(G) is a unitary matrix which transforms the {\g,h*)} basis 
to the {\h* , p, j,t,i)} basis (using the notation defined in the previous section). The discussion above yields 
an explicit decomposition of D(G) into irreducible representations assuming such a decomposition for each 
group algebra C[Z(/i)]. It follows that one can efficiently carry out the quantum Fourier transform over D(G) 
if 

• given h, one can efficiently carry out the QFT over Z{h); and 

• given h and g, one can efficiently express g — tz, where z G Z(h) and t is an element of a fixed 
transversal of Z{h). (Here "fixed" means that may depend on h, but not on g.) 

We remark that the first condition is stronger than merely insisting that there be an efficient quantum Fourier 
transform over the group Z{h): the circuit must be efficiently computable from h in time polylog(|G|). Note 
that the algorithm promised in the second condition implicitly determines a transversal for each Z{h). 
In the following section, we show that these conditions are satisfied for S'„, the symmetric group. 
With such algorithms, the QFT itself is straightforward: 
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1. Convert the basis |.g, h*) to \t, z, h*), where t S Th and z G Z{h). 

2. Rewriting it as \h*,t, z), apply the QFT over Z{h) on the third register conditioned on h* in the first 
register. We obtain the basis \h* ,t,p,i,j), where p S Z{h). Reordering again yields the basis of (3.8) 
above. 

This gives us the QFT over D(G). 
3.3.1 An efficient QFT over D(5'„) 

In light of the discussion above, we show that one can uniformly compute the quantum Fourier transform 
over centralizers in S'„ and, additionally, uniformly factor along a canonical transversal of each centralizcr. 
It follows that we can compute the QFT over D(S'„). 

First, let us determine the centralizer of an element a of 5„ consisting of ci one cycles, C2 two cycles, etc. 
For an element r € Z(a), we have ra = ar and hence rar^^ = a: t fixes a under conjugation. Recall that 
the conjugation action by r is given by the permutation action of r on the labels of the cycle representation 
of a; thus, for each fc, the /c-cycles of a are fixed under the action (ai ... a^) i-> (T(ai) . . . T(afc)). At a 
coarser level, the permutation r fixes, as a set, the elements of 

Wk = {t & {1, . . . , n} \t lies in a fc-cycle of cr} , 

and it follows that Z{a) is the direct product Z{a\) x • • • x .^(ct„), where Cfc is the product of the Ck cycles 
of length k appearing in tr and Z(ak) is the centralizer of Gk (in the subgroup of permutations on Wk\ 
Observe that when (T(a\) . . . T{at)) = (ai . . . at), the element r (restricted to {ai, . . . , at}) lies in the cyclic 
subgroup generated by (ai, . . . , at). It follows that on Wk the element r can be written as the product of 
two permutations, one which cyclicly permutes each fc-cycle, and one which permutes the fc-cycles. This 
subgroup of permutations on Wk is isomorphic to the wreath product Sc^. I T^k — Scf. k Z^*" and we conclude 
that 

n 

Z{a)^^S,,KZl^ . 

Coset factorization Let iJ be a subgroup of G and T a transversal of in G (that is, a set containing 
one representative from each left coset of H in G). {H, G)- coset factorization is the problem of expressing an 
arbitrary element 5 of G as a product th, where h ^ H and t Q T. Observe that if K < H < G is a tower of 
subgroups and Tx and Th are transversals K in H and H in G, respectively, then TkTh is a transversal 
of K in G and (if, G)-coset factorization along this transversal reduces to {K, iJ)-coset factorization and 
{H, G)-coset factorization. 

The induced bases we use to describe the Fourier transform above rely on coset factorization of an element 
along arbitrary, but fixed transversals of the various centralizers Z{a). To see that these can be computed 
effectively in the symmetric groups, let cr G S'„ and 

n 

its centralizer, as discussed above. Our goal is to effectively decompose an arbitrary element tt G 5'„ as a 
product TT = tz, where z G Z[a) and t lies in a transversal T„ (which may be determined implicitly by the 
algorithm). We consider the subgroup chain 

05,, kZ^^ C(i) 05,, k 5^" C(2) 05,,xfe C(3) 5„. 

k k k 

It is an easy exercise to identify natural transversals for the inclusions C(i) and C(3) that permit efficient 
factorization. As for the factorization corresponding to C(2), it suffices to handle individual terms in the 
direct sum of the form 

Sk ^ Si c Ski. ■ 
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For this purpose, we consider an associated action of Ske- We say that a family of sets {Ai, . . . ,Ai) is a 
(fc, ^)-partition if their disjoint union is the set {1, . . . , k£} and each Ai has size k. Then consider the action 
of Ski on the set 

Xk.i = {{Ai,...,Ak} I {Ai,...,Ak) is a (fc, £)-partition of {l,...,k£}} . 
The stabihzer of the element 

a:o = {{!,..., fc},{fc+l,...,2fc},...} 

is precisely a subgroup of the form Si k Sl so we identify a left transversal of Se ix Sl in Ski by identifying, 
for each element x of Xk^i, a permutation that carries xq to x. This can be obtained by "sorting" the sets 
of X according to their smallest element, "sorting" the elements of each individual set, and selecting the 
permutation that carries xq, expressed in the order above, to x, in this sorted order. Factorization along 
this transversal is straightforward by identifying the result of a permutation tt on xq. 

QFTs over the centralizers It is enough to describe how to perform a QFT over wreath products of the 
form W ^Zk wr Si = {IkY x Si. 

The irreps of such groups are easily obtained using Clifford theory (see [10]). Pick an irrep of {I'kY, say 
Lu = (wi, . . . ,uji) where each Wi is an irrep of Zfc, and consider the action of the symmetric group Si on the 
components w^. Define S^j to be the subgroup of Si that fixes this irrep, i.e., the subgroup whose elements 
permute components which are the equal. Thus St is a Young subgroup, equal to a product of symmetric 
groups (each of which acts on a collection of equal indices). If A is an irrep of S^j, then a; (8) A is an irrep of 
= CZkY XI S^. Clifford theory asserts that the irreps of the wreath product are of the form 

t"^^ (w ® A) . 

Woj is traditionally called the inertia group of the irrep uj. 

The QFT over this group can now be constructed using this structure of the irreps. We need to construct 
a transformation which takes us from the basis |(zi, . . . , Zf), tt) d W to the basis |i,a;,A), where t is an 
element of a fixed transversal for in W. In order to do this, we first perform a QFT over {ZkY to 
obtain the basis |w,7r). Now conditioned on \uj), we re-write this basis as \uj,t,iT'), where t is an element of 
a transversal for this Young subgroup and tt' is an element of W^^. Again, conditioned on w, we perform 
a QFT over Wui on the third register to obtain the basis \uj,t,X,i,j), where A is an irrep of S^j and i and 
j label its rows and columns. This step may be carried out by Beals's algorithm [6] for the QFT over the 
symmetric group. This gives the required basis. 

4 Algorithms for approximating link invariants 
4.1 Link invariants 

A knot is a closed (smooth) curve in MP with no self intersections; a link is a finite collection of non- 
intersecting knots. We study knots and links upto deformation; in several cases, we work with oriented 
variants. In particular, recall that a continuous map D : M.^ x [0, 1] —?' M.^ is an isotopy of R'^ if each 
Dt{x) = D{x,t) is a homeomorphism and Da is the identity map. We then identify two knots (or links) if 
there is an isotopy that carries the first onto the second; in this case we say that they are ambient isotopic. 

A link invariant is association of links to algebraic objects (numbers, groups, modules, etc.) that respects 
the equivalence relation of ambient isotopy: equivalent links must be associated with the same object. One 
well-studied framework for studying links is to represent them as certain canonical "closures" of braids and 
explore the algebraic properties of braids that preserve link equivalence in this representation. We consider 
two such closures here: trace closure and plat closure. 
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(a) The trace closure 



(b) The plat closure. 



Figure 2: Braid closures. 



4.1.1 Trace closure 

Trace closure of a braid is obtained by joining the ends (top to bottom) as shown in Figure 2a; Alexander's 
theorem [4] asserts that any link can be obtained in this fashion. Of course, a given link can be represented 
as the closure of many different braids. 

It has been shown by Markov [22] that braids that yield equivalent knots under the trace closure are 
related by the Markov moves. These moves are 



Of] 




rj9, 0,r] € Bn , and 



±1 

ri-li 



A Markov trace is a map 



Bn — > C that is well-behaved with respect to these moves: 
c^iOij) = <f>{r,9) , 



where z = (/)(cr„_i) and z = 4'{cr^]i^). With such a trace map, one can construct a link invariant by defining, 
for any braid 9 that realizes the link, the quantity 



m = (z-z) 



-{n-l)/2 ( t 
Z 



e{e)/2 



(4.1) 



where the linking number, e{9), is the sum of the exponents of the generators Ui appearing in 9. 

It was shown by Tsohantjis and Gould [30] that the quantum doubles of finite groups yield link invariants 
by this approach. Specifically, let (p, h) be an irrep of D(G) acting on the Hilbert space V and let r denote 
the natural representation of the braid group i3„, described in Section 2.2 above, on the space 1/®". Then 
the quantity 

L,^^{9) = {d,.n)-Hh)-^^'^Tr{T{9)) (4.2) 
is a link invariant, where {h) p — Xp{h)/dp and Xp, dp arc the character and dimension of the irrep p of Z{h). 



4.1.2 Plat closure 

A related method to represent links by braids is the plat closure. This is defined on braids with an even 
number of strands where one takes pairs of adjacent strands on the top (and bottom) and joins them as in 
Figure 2b. Similar to trace closure, it can be shown that any link can be represented as the plat closure of 
some braid. 

Birman [7] proved an analogue of Markov's theorem for the plat closure. 
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Theorem 1 (Birman [7]). Given two braids in B2n with the same plat closure, there exists a sequence of 
moves from one to the other such that each move is one of the following: 

• (Double coset move) 9 i — >■ hi9h2, where hi,h2 G H2n (defined below) and, 

• (Stabilization move) 9 i — >■ cr2n9- 

Here, H2n denotes the Hilden group, the subgroup of the braid group i?2n generated by 
(Ti , (J2<j\<J2i and (J2i<y2i-i<J2i+i<J2i ^ Vi e {1, . . . , n - 1} . 

The quantum double of a finite group G likewise yields link invariants by this approach. Consider the 
2n-fold tensor power of an irrep A for which A = A, i.e., A is its own conjugate. For adjacent pairs A®^, 
define the state 

I'i') = 4= E (4-3) 



where the sum runs over all the elements of an orthonormal basis of Va ^nd d^ is its dimension. This 
maximally entangled state is the copy of the trivial representation in the tensor product. Now, define the 
state 

|a) = 1$)®" (4.4) 

and, given a braid 9 G B2n that realizes the desired link under the Plat closure, let the representation of 9 
in the above tensor power be tj>^{9). The plat closure yields a link invariant: 

P1a{9) ^ {a\TA{e)\a) . (4.5) 

Physically, one can view the plat closure as follows: particle-antiparticle pairs are created, their world lines 
are braided and then they are annihilated. The plat closure is the amplitude of obtaining the vacuum 
at the end of this process. One also requires that the particle be its own antiparticle since we need the 
representation to be its own conjugate. 



4.2 Quantum algorithms for link invariants 

In this section, we present quantum algorithms to additively approximate the trace and plat closures of 
braids in some irrep A of D{G). We use the QFT over D(G) to accomplish this: critically, the R matrices 
can be implemented efficiently in the regular, combinatorial representation so we may use the QFT over 
D(G) to implement it inside the irreps. Then we take the trace over the irrep of choice. We describe this in 
more detail below. First, we need the definition of additive approximation (see Bordewich et al. [9]). 

Definition 1. Given any function f : D ^ C and a normalization u : Z+ — > M+, an additive approxi- 
mation for the pair (/, u) is a probabilistic algorithm which given any x ^ D and e > produces an output 
f{x), such that 

¥,[\f{x)-f{x)\>eu{\x\)]<l/^, 

in time polynomial in \x\ and . 

In order to approximate the trace closure of some braid 9 written as a word in the generators at , we first 
construct the representation of 9 in D(G')'^" by implementing the R matrices in D(G). In the case of finite 
group doubles, the R matrix is a unitary operator and, in fact, a permutation operator in the standard basis. 
Recall that it is defined as 

R^J29®9*'^ 
a 

thus its action on a basis element \gih1) (E) \g2h2) is given by 

R\gih*i,g2hl) =^\ggihl,g2{g^^h2)*) = Ih^" gihl,g2h*2). 
g 
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Therefore, the action of tJi = TR is 

0-1 1 31 ^1,52/12) = 132^2:^2' 9iK)- 

This operator is a left multiphcation by a group element followed by a swap which can be efficiently imple- 
mented. 

In order to take the trace, we must be able to produce random basis vectors in the irrep A embedded 
in Vd{g)- The basis vectors inside the irrep {h,p) can be constructed using the fact that it is an induced 
representation. If {wi, . . . ,Wd} are a basis for Vp, then {ti (8) wi, . . .tk ® Wd} is a basis for Vp,/i, where 
ti e Th as before. Assuming that we can embed the vectors \wi) in the vector space Vz(^h)^ then the vectors 
\h*,tj,Wi) are basis vectors of the irrep {h,p) embedded in Vd(g). Now, using the Hadamard test, we can 
find the iimer product 

{v,\FTj,{e)F-^\v,) , 

for any given Vi, where F is the QFT over D(G'). The Hadamard test is a standard tool in quantum computing 
used to estimate the matrix entries of a unitary operator. The description below follows [31] (see also [3]). 

Hadamard test Suppose that J7 is a unitary operator that can be implemented efficiently, i. e., in 0(poly(ri)) 
time and \il>) is a pure state on n qubits which can be prepared efficiently, then one can efficiently sample 
from random variables X, — 1,1}, where 

E[x + zy] = (V'|f/|V')- 



Finally, the Chernoff bound shows that choosing Vi uniformly at random from the orthonormal basis and 
computing the inner product {vi\FTx{0)F~^\vi) is enough to give an additive approximation. 

Chernoff bound If {Xi, . . . , ATfc} are real- valued random variables such that \Xi\ < 1 and E[Xi] = /i for 
all i, then 



Pr 



K 



> e 



<2exp(- — 



In order to apply this to complex numbers Zj = Xj + iYj , we need the union bound along with the Chernoff 
bound. 



Pr 



k 



> e 



< Pr 



k 



> 



V2 



Pr 



k 



> 



V2 



< 4exp 



To determine k (the number of times we need to run the algorithm), let Zi = Xi + lYi where Xi and Yi are 



obtained from the Hadamard test. Let Zi = (dp^h) ^ 



e(e) 



Zi. Then for any e, the probability that the 



average of Zi over k trials differs from the trace closure by more than e is less than 4 exp(— fce^/8). In order 
to make this quantity less than 1/4, we need k > 321og2/e^. 

The algorithm to additively approximate plat closure is similar. We must prepare the state \a) from 
(4.4) efficiently. In order to do this, it is enough to prepare the state |<I>) efficiently. Once we do this, we 
can perform the Hadamard test to estimate the inner product {a\TA{0)\a) and use the Chernoff bound to 
find the additive approximation to plat closure of 9. The state |$) is the maximally entangled state over the 
irrep A(g) A embedded in D(G) (g) D(G). This can be prepared by first creating a maximally entangled state 
over a space of dimension V^^. Then we apply the embedding isometry to embed this state into C[D(G)]®2. 



4.3 Classical algorithm for fluxon irreps of D{G) 

In this subsection, we observe that there is a classical randomized algorithm to additively approximate link 
invariants for certain, combinatorial irreps of D(G). When the strands are colored by irreps of the type 
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A — (/i,tr), where tr is the trivial irrep of Z{h), then it is possible to approximate the link invariants 
classically. In these irreps the R matrix has the action of a permutation matrix. 

R\9i) ® I52) = 1525152^') ® I52) • (4.6) 

This means that for any braid 0, t\{0) is a permutation matrix whose matrix entry can be determined 
efhciently classically (efficient in the number of crossings in 9). Therefore, we can take elements uniformly 
at random from the basis vectors of V^" and compute the inner product with ta{0). Then, by the Chernoff 
bound, we obtain an additive approximation to the trace closure. The difference between the classical and 
quantum algorithms is that, in the quantum case, one has the Hadamard test to efficiently estimate the 
matrix entries of a unitary operator. Here since the R matrix is a permutation matrix, its entries are easy 
to determine classically. Therefore, we have 

Theorem 2. Given a braid b G i?„ of size m, a finite group G and any conjugacy class G of D(G) such that 
there exists a quantum algorithm to additively approximate the trace and plat closures of the braid when the 
strands are colored by the irrep {C,id) ofD{G), in time 0(poly(m,n,\og\G\,l/ e)), where e is the error in 
approximation. 

5 Clebsch-Gordan decomposition 

In this section, we first describe how an efficient algorithm to perform the Clebsch-Gordan decomposition 
over 0(0) for fluxon irreps. Then we give an efficient algorithm for general irreps under certain conditions. 
The following theorem from Curtis and Reiner [10] is useful to understand the Clebsch-Gordan transform. 

Theorem 3. Suppose we have two subgroups H and K of G and two representations p and a of H and K 
respectively. Suppose that we induce these two representations to G, then the tensor product of the two G 
representations can be decomposed into a direct sum of induced representations as 

p \^ ®a t'^^ (p innK-^ ^o- innK-^) \^ , 

d 

where d runs over all {H, K) double coset representatives. 

We can apply this theorem to irreps of D(G) since they are all induced representations from centralizer 
subgroups. Suppose that we have the tensor product of two irreps of D(G), say ([5], p) and {[h],a), then the 
theorem implies that we can write this as 

p t'^ ®(T t'^^ (p iz{g)nz{hr iz{g)nz{hr) • (5.1) 

d 

However, this is still not in the form that we want. In order to obtain it, we need to understand the double 
coset representatives. We show that the {Z{g),Z{h)) double coset representatives also label the different 
conjugacy classes that appear in the product [g] ■ [h] in the conjugacy class algebra. To see this let the 
conjugacy class of g be {g, kigk^^ , . . . , kngk~^} and that of h be {h, lihl^^ , . . . , Imhl^^}, where ki and h 
label the complete set of coset representatives of Z{g) and Z{h) respectively. Now consider all possible 
products of the elements of the two sets. In order to determine the different conjugacy classes that appear in 
the products, we only need to consider elements of the form glihl~^ since anything of the form kjgk~^lihl~^ 

can be conjugated by k~^ to get the former type. 

This is still not enough since two elements ghhl^^ and gljhlj^ could be conjugates. In that case, there 
must be an element z of Z{g) such that z{glihl^^)z^^ = g{zli)h{zli)~^ = gljhlj^. This means that Ij 
is in the same double coset of {Z{g),Z{h)) as k. In other words, the right action of Z(g) on left cosets 
of Z{h) determines the different conjugacy classes which is exactly the different {Z{g),Z(h)) double coset 
representatives. Picking a set of double coset representatives d, we can say that the different conjugacy 
classes that appear in the product [g] ■ [h] are [g/i*^] for all the different double coset representatives d. 
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We now need to determine the number of times each [gh'^] appears in [g] ■ [h] since this determines (some 
of the) multipHcities in the Clebsch-Gordan decomposition. It is enough to determine the number of times 
[gh] appears in the product [g] ■ [h] since its entire conjugacy class would appear the same number of times 
(for any other d the procedure is the same). In order to count this, consider the two groups Z{gh) and 
Z{g) n Z(h). It is easy to see that the latter is a subgroup of the former. If it is a strict subgroup, then 
there exist elements in Z{gh) which do not commute with either 5 or ft, or both. It turns out that in fact, 
any element of Z{gh) either commutes with both g and h or does not commute with both. To see this, 
notice that if s is an element not in Z(g) n Z(h) and it commutes with g but not ft, then we have that 
gh = {ghy — g^h^ — gh^ . But this implies that — h and so s commutes with ft. This means that any 
element in Z{gh) that does not commute with t/, also does not commute with ft and vice versa. 

Now for any non-trivial coset representative s of Z(g) n Z(h) in Z(gh), we have that gh — sghs^^ — g'^h^ ■ 
Since s does not commute with g and ft, we have produced a pair (<;*, ft'*) distinct from {g, ft) such that their 
product is the same. For distinct coset representatives s and t, the pairs (17*, ft*) and {g*,h*) are distinct 
since if not, then g'' = g* and ft'* = ft*. This means that st~^ € Z{g) D Z{h) which is not possible since s 
and t are in distinct cosets of Z{g) fl Z{h). Therefore, for each distinct coset representative s of Z{g) D Z{h) 
in Z{gh), we get a distinct pair {g^, h^) such that g'^h'^ = gh. The same argument holds for any non-trivial 
d. To show that these are all the possible pairs, we make use of the above theorem. When p and a are both 
trivial, Eq. 5.1 takes the form 

([g],tr)®([ft],tr)-0(trtfg;';),(,.)) ■ 

d 

This shows that the number of times [gh"^] appears in [g] ■ [h^] is the index of Z{g) n Z{h^) in Z{gh'^). 

With all this in hand, we can determine the Clebsch-Gordan decomposition. We can re- write Eq 5.1 to 

get 

p t<= ((p iz(g)nz(h)^ izig)nzihr) t''^''''^) , 

d 

by taking the induction in two stages. Now suppose that {p iz{g)nz{h)'' '^'^ iz{g)nz{h)'') t^*'^'* breaks up 
into irrcps of Z{gh'^) as 

{p iz{g)nz{hy'^'^ izigjnzihy) t^''^''''^= ■mp^a.f.p-, 

p.ez(^) 

where mp^a,fi is the multiplicity of the irrep p. We now obtain the Clebsch-Gordan decomposition as 

pt^0at^=0 mp,.,^(/it'=). (5.2) 

fiezi^) 

Since this decomposition is obtained by considering the action of G alone, we need to check if the action of 
ft* is consistent with it. We do this in the next section as we develop the transform. 

5.1 Clebsch-Gordan transform 
5.1.1 Fluxon irreps 

We first describe the Clebsch-Gordan transform for irreps of 0(0) of the type ([(7],tr) (fluxon irreps). For 
this case, we give an efficient transform for any finite group. For two fluxon irreps, the Clebsch-Gordan 
decomposition states that 

([ft],tr)®([5],tr)-0(trt|;^;;\,,))t^-0 m,,,,(/n^). 
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Wc give the transform in two steps, one for each of the above two isomorphisms. In the first step, wc need 
to convert from the basis j/i", g'^) to t), where t is an element of the transversal of Z [g) r\ Z {h''') in 

Z{gh'''). Once we fix the double coset representatives d and the transversal t (for each d), this transformation 
is straightforward. Given g'), we determine their product and the double coset representative d, that this 
product belongs to. Recall that for each element t in T^, we have that g*{h'^Y = gN^ and that each such 
t gives a distinct pair ((/*, (h'^Y). Therefore, we can determine t after ordering the pairs. This gives us the 
first transformation. It is also clear from this that the co-algebra action is consistent across the both sides. 
Indeed, the action of any c* on a state of the form \ti,Vi) (g) \t2, V2) is 



Under the Clebsch-Gordan transform, this state is taken to |f2^i,f3), where Vz is a state determined by Vi 
and V2 ■ This state is given by the second transformation described below. Therefore, the action of c* is 



This shows that the coalgebra action is consistent with this decomposition. 

For the second transformation, we need to block diagonalize the induced representation of the trivial 
irrep of Z{g) n Zih"^) to Z{gh'^). We show how to do this if we can perform eflacient QFT over both groups. 
Suppose that p is an irrep of a group B which is a subgroup of the group A and suppose that we can perform 
a QFT over A and B, then we show how to use these QFTs to give an efficient way to block diagonalize the 
induced representation p t^. We first embed the induced representation into C[^] in the following way. The 
induced representation consists of vectors of the form where t is an element of the transversal of B in 

A and v is a. vector in the irrep space p. We embed \v) into C[i?] by taking an ancilla of size \B\/dp, where 
dp is the dimension of p. This is possible if we know how to perform the QFT over <C[B]. Now that we have 
an embedding into C[j4], we can perform the QFT over A and obtain a basis jp, i, j) where only some irreps 
p appear. Since we know which ones do not appear at all, we can label the irreps of A such that the state 
is |0) \p,i,j) where the first register is of size \H\/dp. We can now discard this register and wc obtain the 
required block diagonalization. Using this technique when p is a trivial representation of i? = Z(yg) n Z{h'^) 
and A = Z[gh'^) gives us a way to decompose the induced representation from the trivial representation of a 
subgroup. This gives us the second transformation and completes the Clebsch-Gordan transform for fluxon 
irreps. Wc have thus shown that if we can perform an efficient QFT over Z{g) and over Z(()) fl Z{h) for all 
g,h ^ G, then we can perform an efficient Clebsch-Gordan transform over D(G). 

5.1.2 General irreps 

For the more general irreps, we show that one can perform efficient Clebsch-Gordan transform if we can 

1. perform QFT and CG transforms over Z{h) and Z{g) n Z{h) for all g,h G G and, 

2. block diagonalize irreps of centralizers restricted to intersections of centralizers. 
The procedure can be split into three steps in the following way. 



c*(|ii,t^i) 1^2,^2)) = ^ K\ti,vi) O h2\t2,V2) = S^t2t^^c}\ti,vi) 1*2,^2) • 



h^hx—c 



C*\t2ti,Vz) = '5{t2ti=c} 1*2^1, -(^s) • 



{[h],p) ® {[g],cr) ^ ((p;z(5)nz(ft<') ®^ iz(g)nz(h'^)) \^^^'''^) f 



.G 



d 




(5.3) 
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In the above, riu is the muhiphcity of v (an irrcp of Z{g) n Z(h)) m the tensor product decomposition of p 
and a restricted to Z{g) n Zihf^) and is the multipUcity of /i (an irrep of Z{gh'^)) when v is induced to 
Z{gh'^) and decomposed into irreps of Z{gh'^). These multipUcities depend on p, a and d in general. 

1. The first step takes us from the basis \h'^,vi,g^,V2) to the basis \{gh'^)'',t, vi,V2), where vi and V2 are 
vectors in the irrep spaces of p and a. This step is the same as in the previous case (for fluxon irreps) 
since we do not operate on the vectors vi and V2. The transversal is picked in the way described above. 

2. The second step can be done if we know how to decompose any irrep of Z(g) and Z{h'^) into irreps 
of Z{g) n Z{hY and then perform Clebsch-Gordan transform over the group Z{g) n Z{hY. This step 
may be done for particular groups. 

3. The third step can be done using the procedure described above for decomposing induced irreps, since 
v is an irrcp of Z[g) n Z{h'^) (the group B above) and it is induced to Z{gh'^) (the group A above). 

This gives an efficient algorithm for the CG transform under certain conditions. In the appendix, we show 
how these conditions are satisfied for D(Zp xi Zg) which is sufficient for universal quantum computation [25]. 

6 Computational complexity of approximating link invariants 

In this section, we show that providing an additive approximation to the plat closure of a braid (when the 
strands are colored by an irrcp of D{G)) is BPP-hard; likewise, we show that providing a multiplicative 
approximation is SBP-hard. For the case of fluxons, however, since we have a randomized classical algorithm 
(presented in Section 4.3) to additively approximate plat closure, this problem is BPP-complete.~ Our results 
also imply that computing the plat closure exactly is #P-complete. In the case of link invariants arising from 
the quantum group Ug(SU(2)) such as the Jones polynomial, it has been shown that approximating them 
additively is BQP-complete and approximating them multiplicatively is ^P-hard. The complexity of these 
1)5(511(2)) invariants appears to be closely related to the fact that the image of the braid group representations 
is dense (in the unitary group). Indeed, Kuperberg [20] shows that denseness is a sufficient condition for 
BQP-completeness, i.e., if the image of the braid group is dense, then this implies the BQP-completeness of 
additive approximations and #:P-completeness of multiplicative approximations. By contrast, as mentioned 
above, the image of the braid group representations arising from D(G) is finite [12], which may explain the 
difference in complexity. A conjecture along these lines is made by Rowell [28]. 

In order to prove our results, we first show that given any randomized computation on n bits, there exists 
a braid on poly(n) strands such that the probability of success of the randomized computation is arbitrarily 
close to the plat closure of the braid. From this, it follows that approximating the plat closure additively 
is BPP-hard and approximating it multiplicatively is SBP-hard. It also follows that computing it exactly 
is =i^P-complete. Additive approximation was defined in Section 4.2. We reproduce it here for convenience 
along with the definition of multiplicative approximation [9] . 

Definition 2. Given any function / : Z) — > C and a normalization u : Z"*" — > R"*", an additive approxi- 
mation for the pair (/, u) is a probabilistic algorithm which given any x Cz D and e > produces an output 
f{x), such that 

T>i[\f{x)~f{x)\>eu{\x\)]<l/A, 

in time polynomial in \x\ and . 

A multiplicative approximation is a special case of a more general value-dependent approximation defined 
by Kuperberg [20]. A probabilistic algorithm to multiplicatively approximate a function is called an FPRAS 
(fully polynomial randomized approximation scheme). 

^When we say BQP-complete (BPP-complete, respectively), we mean PromiseBQP-complete (PromiseBPP-complete, respec- 
tively). There are no known problems which are BQP-complete (BPP-complete, respectively). 



18 



Definition 3. Given / : D — > C, a multiplicative approximation or an FPRAS is a probabilistic algorithm 
which, given any x ^ D and e > 0, produces an output f{x) such that 

PT[\fix)-fix)\>efix)]<l/4, 

in time polynomial in \x\ and . 

If the algorithm is deterministic, then we have an FPTAS (fuUy polynomial time approximation scheme). 
Next, we define the class SBP. 

Definition 4. A language L is in SBP if there exists a polynomial p, a constant e > 0, and a polynomial 
time probabilistic algorithm such that for any input x 

xe L =^ Pr[A{x) accepts] > (1 + e)2-P(l^l) , 

L Pr[A{x) accepts] < (1 - e)2-P(l^l) . 

The class SBP (small bounded probabilistic P) was defined in Bohler et al. [8], which also established 
several fundamental properties of the class. It was shown that SBP contains NP and, in fact, lies between 
MA and AM. From the definitions of multiplicative approximation and SBP, we can see that multiplicatively 
approximating the success probability of an arbitrary randomized computation is (precisely) SBP-hard. 

We now state and prove the main theorem. 

Theorem 4. Let Pl{b) be the plat closure of a braid where the strands are colored by any conjugacy class of 
the alternating group Am, m > 5, which has at least 4 fixed points. Then, additively approximating the plat 
closure is BPP-complete, multiplicative approximation is SBP-complete and exact evaluation is ^P-complete. 

We need the following lemmas. The next lemma is the main result of this section. 

Lemma 5. Given a randomized computation on n bits and any e > 0, there exists a braid on poly{n, i) 
strands and poly {n, ^) crossings with the strands colored by the conjugacy class irrep of the alternating group 
(Am, m > 5) with at least 4 fixed points, such that if the the plat closure of the braid is PI and the probability 
of success of the randomized computation is Ps then we have 

\Ps-Pl\ <e. 

Proof. We use the following model for randomized computation. Let n be the input size. Consider poly(n) 
d-level systems (dits), where d could be exponential in n; we label the levels through c? — 1. Start the 
computation with k = poly(n) dits in the state and i = poly(n) random dits. Perform reversible computa- 
tion on them using the Toffoli gate. We say that the computation is successful if the first dit is a 0. Then 
take one extra dit in the state and XOR the first dit onto it. Now uncompute on the remaining dits to 
get back the initial state and therefore, the computation is successful if the final state is the same as the 
initial state. Let the deterministic reversible circuit be represented by the matrix M and the random dits 
be ri . . . r; . Then the probability of success is 

Pr[Accept] ^ J2 ^^''^i ' ' ' ri\M\0''ri ...re). 

ri,...,ri 

Since M is deterministic, it is a permutation matrix and we can rewrite the probability of success as 
Pr[Accept] = l ^ (OVi . . . r,|M|0'=si . . . s,) = (0|M|,/)) , 

" ri.,....ri 
si,...,se 

where |0) = ^ Er,...r, |0''^i • ■ ■ ^f)- 

In order to relate the success probability to plat closure of a braid where the strands are colored by 
fluxon irreps (or conjugacy classes of G), we use the Ogburn-Preskill encoding [26]. We take as the d level 
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system, a pair of strands closed in a plat, so that the state on the pair is X^gec IS' 9 ^) ~ X^r k) 
the conjugacy class C is of size d. Therefore, pairs of wires closed in a plat can be used as random dits. With 
this encoding, we have that the probability of success is 

Pr[Accept] = ^ {ici,c^^)''9i9T^ ■ ■ ■ 9e97^\M\{ci,c^^)''hih^K . . hihj^} , 

where the g's and the /I's run over the entire conjugacy class and |c, c~^) represents the |0) state. These can 
be thought of as variables and the cZits in the zero state as group constants since they need to be in a specific 
state. If there were no need for dits in the zero state and if M can be realized as a braid, then it is easy to 
see that the above expression is exactly the plat closure. In order to realize M as a braid we need to have 
access to group constants. Therefore, group constants are needed for the circuit as well as for the starting 
state. Now, if group constants are available, then M can be realized as a braid if the group is sufficiently 
strong (such as a simple group). It has been shown by Maurer and Rhodes [23] that in simple groups, any 
function can be realized as a word in the variables and group constants. To construct only the Toffoli gate, 
less rich groups might work as well (see [5]). 

In order to produce group constants with plat closure, our strategy is to introduce more strands and 
braid them so that when they are closed in a plat, the only solution to the resultant set of equations over 
the group, is the set of group constants we need. We can then normalize the plat closure so that this part 
of the braid contributes unity. Assume that the constants needed are ci, c^^ , . . . Cm, , (ci, cj~^)'^, where 
ci, . . . , Cm generate the group (we assume that G has a conjugacy class which contains its own inverses and 
can generate the group). We now show how to generate constants. 

The invariant at a fluxon irrep Recall the Wirtinger presentation of the knot group IC{K) — Tri{R'^\K) 
of a knot K (or link): beginning with a knot diagram** for K, one introduces an orientation for the knot and 
a generator Xs for each stroke of the diagram (that is, each connected arc of the diagram, such as the red 
component of the knot diagram of Figure 3b). With each crossing of the diagram, such as the one pictured 
in Figure 3a involving the labeled strokes Xt, Xa, and x^, one introduces the relation 

Xb = Xi^XaXt. (6.1) 

(One can recover this relation by associating each variable Xs with the path from a fixed base point that 
travels around the stroke s in a direction consistent with the right-hand rule and recording the effect of 
"pushing" a loop underneath the top strand.) The quotient of the free group generated by the {xg] by these 
generators is isomorphic to the knot group. 

As the knot group is an invariant of the knot, one can immediately obtain other (weaker) invariants as 
functions of the knot group. We focus on 

h{K, G) = \H{K, G) I where H{K, G) = {0 : tC{K) ^ G | </> a homomorphism} , 

the number of homomorphisms of K,{K) into the group G. With the Wirtinger presentation describved above, 
the quantity h{K, G) can be expressed as the number of maps h : {xs\ G that satisfy the relations (C.l) 
in the sense that 

(l>[xc) = (/)(xt)"V(a;c)0(a;t) 

for each crossing. Observe that any such map necessarily carries the generators {xs\ into a single conjugacy 
class of G, and we refine the notation above by focusing on those maps associated with a particular conjugacy 
class G of G: 

h{K; G, G) = \H{K- G, G)| where H{K- G, G) = {</> : 1C{K) ^ G | a homomorphism, Vx^, 0(x,) e G} , 

^Recall that a knot diagram is a 2-dimensional projection of a knot in general position, so that no three lines intersect at a 
point, with explicit annotations that determine, for each crossing, which line occludes the other. Figure 3b is a knot diagram 
of the unknot. Figure 3a is an example of the typical indication of a crossing. 
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which we abbreviate h{K, C) (and H{K, C)) when G can be inferred from context. 

As discussed above, a braid b G i?2n induces a Hnk L{b) via the plat closure, and any (efficiently presented) 
knot can be given such a description (as the knot induced from the plat closure of a braid) in a straightforward 
fashion. For a group G, and the fluxon representation A = (C, 1) of D{G), this yields the knot invariant 
P1a{0) (where is a braid that yields the knot under the plat closure). For this combinatorial case (arising 
for a fluxon representation), we have 

PlA = hilCiK);G,G). 
With this equality in place, we shall argue about h{IC{K); C) rather than Pl\. 



(a) A crossing in a (oriented) knot diagram. (b) A stroke in a knot diagram. 

Figure 3: Knot diagrams, strokes, and the Wirtinger presentation. 



Generating groups constants As mentioned above, we assume that G is a finite group and that C — 
{h^ I g G G} is a conjugacy class that generates G. We begin by showing how to construct a knot that — with 
controlled error — permits us to distinguish a particular set of generators of G (up to automorphism) . 
Let c = (ci, . . . , Ct) denote a sequence of elements of C that generate G and 

[c] = {0(c) = (0(ci), . . . , ^{ck)) I e Aut(G)} . 

Our goal is to construct a knot kc, together with a distinguished collection of strokes (si,...,Sfc), with 
the property that uniform selection of a homomorphism -ijj from H{k,C), with overwhelming probability, 
yields one for which tp{s) = (V'(si), . . . ,^/'(sfe)) € [c]. (Observe that Aut(G) acts on the collection of legal 
homomorphisms, so our demand — that the resulting knot distinguish a particular orbit under the Aut(G) 
action — is the strongest requirement of this type on which we can insist upon.) 

The construction begins with a collection of k (oriented) circles; see Figure 5a. Without further con- 
straints, the knot group of this link is free: legal maps 4' may assign arbitrary elements of C to each circle. 
We introduce constraints among the circles with band connections; see Figure 5b. A band is a pair of strands 
that are constrained to be equal (under any homomorphism into G), but are oppositely oriented. They are 
convenient because of the simple effect they have upon crossing other strokes: 

1. Should a band cross beneath another stroke, as in Figure 4b, the constraints of (3a) preserve the band 
condition: if the pair of strokes on the left side form a band (are constrained to take equal values under 
any homomorphism), the same is true of the two strokes on the right. Furthermore, the values assumed 
by these two bands are related by conjugation of the value of the stroke under which they cross. 

2. Should a band cross over another stroke, effectively dividing it into three strokes as in Figure 4b, the 
constraints of (3a) demand that the "incoming" and "outgoing" values assumed by the stroke so crossed 
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over are the same, as though the crossing had never happened. The smah interior stroke will take on 
another (conjugate) value, but we shall never allow these to interact with other strokes. 



y 




y 






















X 




— ^ 



yxy 



(a) A band. 



(b) A band crosses beneath a 
stroke. 

Figure 4: Bands and band crossings. 



(c) A band crosses over a stroke. 



Since ci, . . . , c„i are conjugates, we can write each Ci as g ^cig, for a group element g. As the Ci generate 
the group, however, g may be written as a word in the Ci and we have, for each i an equation 

c^ = , (6.2) 

where Vi is a word in the Ci and their inverses. Introducing a family of indeterminates xi , . . . , and replacing 
each appearance of Ci in (6.2) with the variable x^, we conclude that the q satisfy the equations 

Xi = Xi^ , (6.3) 

where Wi is the word in the Xi and their inverses obtained from Vi . We now adopt a construction of Johnson 
[17], who shows how to construct a knot that imposes these relations. Each such equation Xi — x]"* is 
imposed by introducing a band between the circle associated with xi and the circle associated with Xi (see 
Figure 5h) and passing it through the circles associated with the variables appearing in Wi ; by passing such 
a band "through" such a circle — one over and once under as in Figure 5c — the generator associated with 
the circle operates on the band value by conjugation while leaving the generator associated with the circle 
unaffected. In general, the words Wi may contain references to the variable Xi, in which case one passes the 
band through the circle associated with Xi, as in Figure 5d. 

This construction produces a knot we call K^. By construction, associating the strokes the k original 
circles with the variables Xi, there is a homomorphism -0 of IC{K^) into C <Z G for which tjj : Xi ^ Ci. As 
remarked above, if € Aut(G) which fixes the conjugacy class C, then o i/; is also an element of H{K^, C). 
However, there will, in general, be other solutions to the equations (6.3): in particular, the map which carries 
all strokes of the knot to a particular, fixed element c € C satisfies the equations (6.3). Our goal is to add 
some strands to the knot so as to "amplify" the solution of interest — the map 4^ ■ 

X s I y Ci so that it appears 

with overwhelming majority if an element of selected at random from H{-,C). 

To minimize the contribution from solutions not related by an automorphism, consider a undesirable 
solution d — (c?i, . . .,dk) (for which there is an homomorphism ip' : Xi i— >■ di and yet d ^ [c]). We now 
introduce a new circle, associated with the variable ynow, and as before introduce an equation of the form 

y,,„ = x-^^^y- , (6.4) 

where 'w{x) is a word in the Xi. In order to "suppress" the homomorphism associated with d, we shall choose 
w in such a way that 
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(a) Initial circles. (b) A band between two circles. 




(c) A simple relation: = y. (d) The equation a;^^ = y. 



Figure 5: Building a set of generators. 

1. w{d) = 1: thus any homomorphism -0' which carries Xi to di must necessarily carry t/ncw to the same 
value as xi = di. (Specifically, the only solution to the equation j/ncw = x™'''^''*'"™ = 2:5'°°" is j/ncw = 2^1-) 

2. w{c) = a, a group constant for which the equation j/now = has multiple solutions, fn this case, 
the homomorphism tp : Xi ^ Ci may be extended in multiple ways with an assignment to j/ncw (any 
solution to 2/ncw = a;"*'"™ will suffice). 

Lemma 6 below guarantees that there is an element a G An,n > 5 satisfying item 2 above. Lemma 7 
below guarantees that when G is simple (i.e., has no nontrivial normal subgroups), there is a word w in the 
variables xi, . . . ,Xk so that w{d) = I and w{c) — a for any fixed a € G. Focusing now on a specific group 
G — An, n > 5, let a be the element which maximizes the number of solutions to the equation y = c"". It 
follows immediately that, for the new knot K' induced from this process, 

|{0 e H{IC{K');C, G) I ^{x,) = q}| >2\{^e H{IC{K');C, G) \ i^{x,) = d,}\ . 

Repeating this process i times for each of the offending vectors d induces a knot so that for any d ^ [c] 
we have 

\{<P e H{]C{Ki)-C, G) I ^{x^) = q}| > 2^ |{0 e H{]C{Ki)-C, G) ^ G \ i^{x,) = d^}\ . 
(Observe that, as presented, this is only efficient if G = An has constant size.) 

Lemma 6. Consider the alternating group on m symbols A„i and consider any conjugacy class with at 
least 4 fixed points, say, m — 3, m — 2, m— 1 and m. Since the conjugacy class is non-trivial, it has a 
k cycle with k > 2. Pick a in this class to have (12 ... fc) in that k cycle. Now, consider the equation 
y = Q!^^ This equation has at least two solutions. 

Proof. It can easily be checked that two solutions are given by elements that differ from a only in the k 
cycle. The first one contains (12 ... k — Im) and the second (1 2 ... A; — 1 m — 1) in that k cycle. □ 
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Lemma 7. Let G be a finite group and Fk be the free group on k generators {xi, . . . , Xk}- For an element 
d = (di, . . . ,dk) G G'' , let (pd be the evaluation homomorphism that carries Xi to di and define 

Ad = ker^d = {w e Ffc | = 1} . 

For two elements c,d£ G^ , define G to he the subgroup generated by the Ci, D to be the subgroup generated 
by the di, and 

Ad{c) = {Mw) \ weAd}. 

Then 

1. Ad{c) is normal in G. In particular, when the Ci generate G, Ad{c) is normal in G. 

2. If c and d are not related by a homomorphism of G ( which is to say that no homomorphism from D to 
G carries each Ci to di) then Ad{c) ^ 1. 

Proof. For part 1, ker((/)(;) is certainly normal in Fk since (pd is a group homomorphism. Let C be the 
subgroup of G generated by the and g e C; we wish to show that g^^Ad{c)g — Ad{c). If w e F\^ is a 
word in the free group for which if)c{w) = g, we have w~^AdW = Ad and hence that 

g-^Adic)g = (l)c{w-^Adw) = M^d) = A<j(c) , 

as desired. 

As for part 2, let C denote the subgroup generated by {ci} and D the subgroup generated by {di}. 
Observe that if Ad{c) is trivial we have ker^^ C ker(/)c- In this case, the natural quotient map 

q,:D^F/keT^d^C = F/kei(j), 

yields a homomorphism -0 : I? — > C for which ip : di ^ Ci. To be precise, let Qc denote the quotient map 
qc '. F/hevcjid — ^ F/hevipc and let id denote the inverse of the isomorphism (j)d induces from F/kei(j)d to D; 
observe that id{di) = Xi (ker (pd)- Then the map ip — (j)c° Qc^ id a homomorphism of D onto G that carries 
each di onto c^; see Figure 6. 



^ G 




Figure 6: The homomorphism tp : di Ci. 

□ 

The above lemma implies that if G is simple then Hi — G {Hi cannot be trivial because c and d are not 
related by an automorphism). This means that there exists a word w in variables Xi such that w{c) = z and 
w{d) = 1. □ 

Proof of Theorem (4). The first two follow from the definitions of BPP and SBP and the following lemmas. 
The #P-completeness follows from the fact that exact evaluation of the success probability of a randomized 
computation is #P-complete. □ 
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7 Conclusions 



In this paper, we gave an efficient circuit for the quantum Fourier transform over D{G), the quantum double 
of a finite group G and show how to apply it to D(5'„). We used this circuit to give efficient algorithms for 
approximating link invariants. We then gave an efficient circuit for the Clebsch-Gordan transform for fluxon 
irreps of D(G) and, under certain conditions, for general irreps. We gave an example of a quantum group, 
namely D(Zp xi Zg) (which is powerful enough to do universal quantum computation) for which we show how 
to perform the QFT and CG transforms. We also showed how to simulate topological quantum computation 
inside exponentially large irreps of D(G) efficiently using the Fourier and Clebsch-Gordan transforms. We 
then showed some hardness results for approximating and exact evaluation of link invariants arising from 
D(G). We showed that additive approximations of link invariants arising from irreps of D(G) are BPP- 
hard and multiplicative approximations are SBP-hard and exact evaluations are #P-hard. We also gave an 
efficient randomized algorithm to additively approximate the link invariants when the conjugacy class (or 
fluxon) irrep is used. This shows that for the fluxon irrep the problem is BPP-complete. 

In our hardness results, we need that the size of G be constant. An immediate question is how to extend 
this to asymptotically growing group sizes. One way is to make the procedure to kill unwanted solutions in 
our proof more efficient. Another question is - to which groups can these hardness results be extended. It 
is possible that the hardness results proved here for An are also true for certain non-solvable groups which 
are not simple. For example, it is immediately true for Sn by restricting to conjugacy classes which are also 
in An- But it could be true for a more general class of non-solvable groups. Finally, it would be interesting 
to determine the power of a D(G) computer where we do not use post-selection (with post-selection, it is 
universal for quantum computation). Another interesting question regarding post selection is, what class of 
operations or gates can be extended to universal quantum computation when post selection is used. 
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A Fourier and Clebsch-Gordan transforms over D{Zp xi Zg) 

It is shown in [25], that the group G ~ Zp » Zq can be used to perform universal quantum computation 
when p and q are prime and q\(j) — 1). Here we show how to perform the Fourier and Clebsch-Gordan 
transforms over D(G). First, in order to fully describe the group, we have to pick the homomorphism from 
Zg to Aut(Zp). Such homomorphisms are characterized by elements a such that a'^ = mod p. Having 
picked such an a, we see that the group multiplication in G is (ai, 6i)(a2, 62) — (11 + a2a^^ , bi + 62)- Notice 
that when a = 1, we obtain the direct product of Zp and Zg. 

We now describe its irreps. There are q one dimensional irreps and (p ~ i)/q, q dimensional irreps. The 
one dimensional irreps are obtained as the extension of the trivial irrep of Zp to G and tensored with each of 
the q irreps of Zg. The q higher dimensional irreps are obtained as induced representations from a non-trivial 
irrep of Zp. The irreps of Zp are characterized by k and are of the form exp(27rifca/p) = Wp" for a group 
element a. All those k in the orbit of ka^ for b £ Zq will be induced to the same irrep of G. The induction 
can be given as 

Pfc(a,5)= ^c.p^'^"">)(s-6|. 

We now describe the centralizers of G. The centralizer of the identity element (0, 0) is G. The centralizer 
of any element (a, 6), where a 7^ is Zp. The centralizer of (0, 6), where 6 7^ is Zg. 

The Fourier transform is now easy to construct. Notice that since two of the centralizers are abelian, 
their QFT is efficient. For the QFT over G, first take the group basis \a,b) and convert it into the basis 
\z, t), where z € Zp and t is a transversal which can be picked to be an element of Zg. Now perform a QFT 
over Zp in the first register. Then conditioned on the value in the first register perform a second transform 
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as follows. If the value in the first register is a non-trivial irrep of Zp, then do nothing since we already 
have an irrep of G. If the first register has a trivial irrep of Zp, then perform a QFT over Zg in the second 
register. This gives us a QFT over G. Since we can now perform efficient QFTs over all the centralizers, we 
can perform efficient QFT over 0(0). 

For the Clebsch-Gordan transform, as the first condition, we need to perform CG transforms over all 
centralizers and their intersections. Since all intersections and two of the centralizers are abelian, their CG 
transforms are efficient. For the CG transform over G, we only have to consider the case of high dimensional 
irreps. For these irreps, the tensor product looks like 

s,t 

The Clebsch-Gordan transform, if fc -I- ^ 7^ 0, is 

|s,t) \t-s,t). 

This produces the state 

s.t s,t 

which is a direct sum of q irreps pka'+i for s G Zg. Now, if fc + / = 0, then again we first perform the above 
transform to get 

s,t 

Notice that when s = 0, this is the regular representation of Zg. Therefore, we have to perform a conditional 
QFT (conditioned on s = 0) additionally. Then we would obtain a direct sum of the q one dimensional irreps 
and q — 1 high dimensional irreps Pk(a'-i) for s G Zg, s ^ 0. This completes the CG transform over G. We 
can now perform QFT and CG transforms over centralizers and intersections of centralizers. 

To complete the conditions for CG transform over D(G), we need to block diagonalize restricted to 
intersections of centralizers i.e., block diagonalize pk when restricted to Zp and Zg. From the structure of 
Pk, it is easy to see that when restricted to Zp, it is already diagonal and when restricted to Zg, it is the 
regular representation of Zg. Therefore, a QFT over Zg would diagonalize it. Thus we can perform efficient 
QFT and CG transforms over D(G). 

B Quantum computation with Anyons 

Anyons are particles which exist in two dimensions and have exotic statistics. Anyons are useful for quantum 
computation because quantum information can be stored on a system of anyons in a non-local fashion. This 
means that local errors do not corrupt the quantum information and a computer based on anyons will be 
inherently fault tolerant to local errors. For a tutorial on quantum computation using anyons, see the notes 
by Preskill [27]. In this section, we give short introduction to anyons described by D(G) for a finite group 
G. Then we show how one can use the QFT over D(G) to simulate an anyon computer efficiently. This has 
been shown in [27, 24], but it was assumed that the dimension of the irreps of D(G) are of constant size. 
Here we use the QFT over D(G) to simulate an anyon computer in potentially large irreps of D(G). 

The Hilbert space of an anyon (whose symmetries are described by D(G)) is an irreducible representation 
of D(G). Recall that irreps of D(G) are characterized as {h,p), where ft. is a representative element of a 
conjugacy class and p is an irrep of Z{h). The group elements h are called fluxes and the irreps p are called 
charges. An anyon, in general, has both flux and charge. Anyons which transform as {h,ti), where tr is the 
trivial irrep of Z(h), are called fluxons. On the other hand, if we pick h to be the identity element of the 
group, then Z{h) = Z{e) = G. Anyons described by the irreps (e,p), where p is an irrep of G, are called 
chargeons. Recall from (3.9) that the action of the R matrix on a pair of anyons is 

R\9i,vi) «> 152,^2) = \g29ig2^,p{k^^\_^ ig2kgjvi) (g) 152,^2)- (B.l) 
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In the special case of fluxons, this action reduces to 



R\gi)®\92) = \g2gi92^)®\g2)- (B.2) 

This is the action when we wind the first anyon around the second in the anticlockwise direction. The braid 
operator TR is given by 

B\gi) ® \g2) = 132) ® |<72gi.92"')- (B.3) 

The action on chargeons can be determined similarly. Given an irrep {g,p), the conjugate irrep is {g~^,p), 
where p is the conjugate irrep of p. The action of p is simply the complex conjugation of the action of p. 
Given two anyons in conjugate irreps, the state of trivial total flux and charge is the maximally entangled 
state 

Vl^l"P g,v 

where the sum is over all g in the conjugacy class C and all v in the irrep p. 



Simulation of anyons 

As described in [27], in order to perform universal quantum computation, we need to be able to 

1. Prepare any state in the Hilbert space of a pair of anyons which correspond to conjugate irreps. 

2. Perform braiding of anyons around each other and around ancillas. 

3. Fuse pairs of anyons and measure the flux and charge of the resulting particle. 

We now show how one can simulate each of these steps using the Fourier transform and the Glebsch-Gordan 
transform over D{G). Initial state preparation can be carried out if we can construct states which lie inside 
pairs of irreps of D{G). To do this, we can use the Clebsch-Gordan transform. First we embed the state 
in the direct sum of the CG decomposition and then perform the inverse CG transform. Next, in order to 
perform braiding on the states, we need to implement the R matrix inside the irrep. We use the same trick 
as before and implement the R matrix in the regular representation of D{G) and use the QFT. Finally, in 
order to simulate fusion, we again make use of the GG transform. 
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